Single-particle dissipation in TDHF studied from a phase-space perspective 



N. Loebl 1 , A. S. Umar 2 , J. A. Maruhn 1 , P.-G. Reinhard 3 , P.D. Stevenson 4 , and V. E., Oberacker 2 
1 Institut fuer Theoretische Physik, Universitaet Frankfurt, D-60438 Frankfurt, Germany 
2 Department of Physics and Astronomy, Vanderbilt University, Nashville, Tennessee 37235, USA 
3 Institut fuer Theoretische Physik II, Universitaet Erlangen-Nuernberg, D-91058 Erlangen, Germany and 
4 Department of Physics, University of Surrey, Guildford, GU2 7XH, UK 
(Dated: August 30, 2012) 

We study dissipation and relaxation processes within the time-dependent Hartree-Fock approach 
using the Wigner distribution function. On the technical side we present a geometrically unrestricted 
framework which allows us to calculate the full six-dimensional Wigner distribution function. With 
the removal of geometrical constraints, we are now able to extend our previous phase-space anal- 
ysis of heavy-ion collisions in the reaction plane to unrestricted mean-field simulations of nuclear 
matter on a three-dimensional Cartesian lattice. From the physical point of view we provide a 
quantitative analysis on the stopping power in TDHF. This is linked to the effect of transparency. 
For the medium-heavy 40 Ca+ 40 Ca system we examine the impact of different parametrizations of 
the Skyrme force, energy-dependence, and the significance of extra time-odd terms in the Skyrme 
functional. 

PACS numbers: 21.60.-n,21.60.Jz 



I. INTRODUCTION 

Time-dependent Hartree-Fock (TDHF) theory pro- 
vides a fully self-consistent mean-field approach to nu- 
clear dynamics. First employed in the late 1970's [IH1] 
the applicability of TDHF was constrained by the lim- 
ited computational power. Therefore, early applications 
treated the problem in only one spatial dimension, uti- 
lizing a very simplified parametrization of the nuclear 
interaction. Due to the increase in computational power, 
state-of-the-art TDHF calculations are now feasible in 
three-dimensional coordinate space, without any symme- 
try restrictions and using the full Skyrme interaction jSJ- 

HQ]. 

In this work the Wigner distribution function |11) is 
calculated as an analysis tool to probe the phase space be- 
havior in TDHF evolution of nuclear dynamics. In com- 
parison to previous work |12) , where the Wigner analysis 
was performed in one and two dimensions, we are now 
able to carry out both the TDHF simulation and the 
phase-space analysis in three dimensions. Transforma- 
tion from coordinate-space representation to phase-space 
representation, i.e. calculating the Wigner distribution 
from the density matrix, still remains a computation- 
ally challenging problem. Here, we present a fully three- 
dimensional analysis which allows the study of relaxation 
processes simultaneously in all directions in fc-space. An 
early one-dimensional study of the Wigner function for 
TDHF can be found in [T5] . 

The paper is outlined as follows: In Sec. |TT] we intro- 
duce the Wigner distribution function and discuss the 
numerical framework used in this work. We then intro- 
duce the principal observables summarizing the local or 
global momentum-space properties of the Wigner func- 
tion. First the quadrupole operator in momentum space 
which gives rise to the usual deformation parameters j3 
and 7 to probe relaxation processes in dynamical calcu- 



lations. In addition, we define an estimate for the occu- 
pied phase-space volume to obtain a relation between the 
fragment separation in momentum and coordinate space. 

Sec. |III| provides a detailed discussion of the central 
40 Ca+ 4U Ca collision, paying particular attention to the 
effect of transparency. We discuss the impact of different 
Skyrme parametrizations on the relaxation behavior, as 
well as the dependence on the center-of-mass energy for a 
fixed Skyrme interaction. We also examine the influence 
of extra time-odd terms in the Skyrme functional. 



II. OUTLINE OF FORMALISM 

A. Solution of the TDHF equations 

The TDHF equations are solved on a three- 
dimensional Cartesian lattice with a typical mesh spac- 
ing of 1 fm. The initial setup of a dynamic calculation 
needs a static Hartree-Fock run, whereby the station- 
ary ground states of the two fragments are computed 
with the damped-gradient iteration algorithm [14l US] . 
The TDHF runs are initialized with energies above the 
Coulomb barrier at some large but finite separation. The 
two ions are boosted with velocities obtained by assum- 
ing that the two nuclei arrive at this initial separation 
on a Coulomb trajectory. The time propagation is man- 
aged by utilizing a Taylor-series expansion of the time- 
evolution operator |16j up to sixth order with a time step 
of t = 0.2 fm/c. The spatial derivatives are calculated us- 
ing the fast Fourier transforms (FFT). 



B. Computing the Wigner function 

The Wigner distribution function is obtained by a par- 
tial Fourier transform of the density matrix p(r— | , r+-§ , t) , 
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with respect to the relative coordinate s = r — r' 
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Because fw is not positive definite, it is misleading to 
consider the Wigner function as a phase-space probability 
distribution. We will refer to the appearance of negative 
values for /w in Sec. [TTT 



Evaluating the Wigner function in six-dimensional 
phase space is still a computational challenge and only 
possible employing full Open MP parallelization and ex- 
tensive use of FFT's. The determing factor is the grid 
size, which results in 

N 2 log (N x ) * N 2 log (N y ) * N 2 log (JV, ) (3) 

steps to provide the Wigner transform in full space, where 
N x ,N y , N z are the grid points in each direction. Stor- 
ing the Wigner function reduced to the reaction plane, 

i- e - /w \ x iV = 0, z, k) at one time step will consume 
~ 140 Mb of disk space for the applications presented in 
Sec. |III| Going to larger grid sizes, needed for heavier sys- 
tems, and/or storing the full three-dimensional Wigner 
function will clearly result in entering the Gb regime. 



C. Observables 

In this section we discuss some of the observables used 
in our analysis. In order to avoid any misunderstand- 
ings we will label all observables evaluated in momentum 
space with a subscript k, and all observables in coordi- 
nate space with a subscript r. 



1. Quadrupole in momentum space 

As an observable to probe relaxation in phase-space 
quantitatively, we evaluate the quadrupole operator in 
momentum space. The local deviation of the momentum 
distribution from a spherical shape is a direct measure 
for equilibration. The local quadrupole tensor in fc-space 
is given by 

Q£(r,f)= / d 3 fc[3(fc. i (r,0)(fc,(r,t))-(k 2 (r,t)) ( 5 u ] , 



. (4) 

using the m-th moment from the local momentum distri- 
bution 

(v(m) , _ Jd3fc(k-(k(r,t))r/^(r,k,t) 

with (k(r,i)) denoting the average local flow 
Jd 3 fck/^(r,k,i) 



<k(r,t)> 
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(6) 



The spherical quadrupole moments Q 20 (r,t) and 
Qf 2 (r,t) are computed by diagonalization of Q % £(r,t) 
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with A3 > A2 > Ai labeling the eigenvalues of Q l ^(r,t). 
Switching to polar notation the observables 



&(r,t) = /&M + 2^M 

V2(3 22 (r,t) 180°, 

7fc(r,t) = arctan— — — — 

P2o(r,i) 7T 

are obtained via the dimensionless quantities 
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where 



r k (r,t) = V(k(r,i)) 2 Mr,t) 



(9) 
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accounts for the local rms-radius in fc-space. The norm 
is defined such that 

p(r,t) = J d 3 fc/^(r,k,t). (14) 

In the presented formalism it is straightforward to de- 
fine global observables. The global quadrupole tensor is 
calculated by spatial integration 



J d 3 rp(r,t)Q^(r,t). (15) 



Applying the same diagonalization as in the local case ([7]) 
we end up with a global definition for 0fP{t) and /3f 2 (i). 
For the following results we exclusively use the global 
definition since it is more compact and allows the simul- 
taneous visualization of multiple time-dependent observ- 
ables. 



2. Quadrupole in coordinate space 

To illustrate the global development of a reaction, we 
will also use the expectation value Q 2 . = (Q 20 ) of the 
quadrupole operator in coordinate space. 



3. Occupied phase space volume 

To give a rough measure for the phase-space volume 
occupied by the fragments during a heavy-ion collision 



160 MeV 



we assume a spherical shape of the local momentum dis- 
tribution. Adding up the fc-spheres 



^(r,t) = ^(k 2 (r,i)) 3 / 2 



(16) 



leads to the total occupied phase-space volume 

V k (t) = J d 3 rV k (r,t). (17) 

III. RESULTS AND DISCUSSION 




It is the aim of this work to provide a quantitative 
analysis of the magnitude of relaxation processes occur- 
ring in TDHF. Therefore we will vary a single reaction 
parameter, while all the other parameters are fixed. The 
40 Ca+ 40 Ca-system provides a suitable test case. This 
particular choice is motivated by Ref. [T7] where the ap- 
plicability of TDHF was demonstrated up to very high 
energies. All calculations in this section were done for 
central collisions (impact parameter 6 = 0). The numer- 
ical grid was set up with 36 x 24 2 grid points. 



A. Variation of the Skyrme force 

In the first set of calculations we vary the Skyrme 
parametrization. Figure [T] shows the results of a central 
40 Ca+ 40 Ca collision with the Skyrme parametrizations 
SLy4, SLy6 [E], SkMs pj5], SkI3, and SkI4 [20]. While 
SkMs was chosen as an example for an outdated interac- 
tion, the SLy(X) set of forces was originally developed to 
study isotopic trends in neutron rich nuclei and neutron 
matter with applications in astrophysics. The Skl(X) 
forces take the freedom of an isovector spin-orbit force 
into account. This results in an improved description of 
isotopic shifts of r.m.s. radii in neutron-rich Pb isotopes. 

The global development of the reaction is visualized 
in subplot (d). The time-dependent expectation value 
Q 20 (t) shows the five trajectories initially in good agree- 
ment but finally fanning out. A similar splitting behavior 
depending on the employed Skyrme parametrization was 
already found in [31]. While the two Sk(X)-forces show 
a full separation of the two fragments, there is a slight 
remaining contact between the fragments for the case of 
SLy6, which will result in complete separation in a longer 
calculation. However the trajectories for SLy4 and SkMs 
show a merged system in the final state, which was found 
to persist in long-time simulations. 

We now consider the relationship between the observed 
characteristics in coordinate space with the dynamics in 
phase space. Subplot (a) shows the /^-value, measur- 
ing the global deviation of the momentum distribution 
from a sphere. The initial /3fc-peak is strongly damped 
for all five Skyrme-forces. While the time development 
for all parametrizations remain in phase up to the second 
peak, later it starts to vary and continue with damped 
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FIG. 1. Global observables (3 k (t) (a), y k (t) (b), V k (t) (c), 
and Qf{t) (d) are shown for a central 40 Ca+ 40 Ca collision 
with a center-of-mass energy -B c . m . = 160 MeV. Each curve 
corresponds to a different Skyrme force as indicated in the 
legend. 
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TABLE I. Effective masses of Skyrme parametrizations used 
in this work are listed in connection with the maximal /3k- 
values from the plots (a, e) in Fig. [I] 



Skyrme force 


m* /m 




SkM 


0.79 


0.0116 


SLy4 


0.70 


0.0111 


SLy6 


0.69 


0.0106 


SkI4 


0.65 


0.0102 


SkI3 


0.57 


0.0095 



oscillations. For a better visualization the first peak is 
magnified in subplot (e). The taller the /3fc-peak the 
longer the fragments will stick together in coordinate 
space. The effect appears to depend on the effective mass 
to* /to. Smaller effective masses give rise to a smaller /3^- 
deformation. Table [I] summarizes the m*/m- values as- 
sociated with the maximal deformations /3™ aa: for all the 
Skyrme forces used in this work, ft seems logical that the 
to* /m-dependence is visible in the phase-space analysis 
since it is directly linked to the nucleons' kinetic motion, 
ft is harder to randomize the directed motion of a nu- 
cleon with a higher effective mass than it is for a nucleon 
with a smaller TO*/TO-dependence. 

Subplot (b) shows the 7fc-value which indicates, 
whether a deformation is prolate, oblate, or triaxial |22j . 
For the present scenario of a central collision the 7fc-value 
jumps between prolate and oblate configurations indicat- 
ing that the momentum distribution oscillates between 
being aligned primarily in the beam direction or trans- 
verse to it. For the sake of completeness we additionally 
present the occupied phase-space volume (c) which will 
prove more useful for the next reaction parameter to be 
discussed: the center-of-mass energy. 



B. Variation with the center-of-mass energy 

As a second reaction parameter the center-of-mass en- 
ergy, E c m , is varied. Results are presented for ener- 
gies ranging from E c , m , = 2 MeV/nucleon up to -Ec.m. = 
3 MeV/nucleon. The Skyrme interaction now is fixed to 
be SkM. For the case of the lowest (highest) energy 
Video [I] (Video [2]) provides a video visualizing the re- 
action in phase space. The calculation done with the 
lowest energy £ c . m . = f 60 MeV shows two fully sepa- 
rated fragments in the exit channel. In contrast, the case 
with the highest energy (as well as the one at an inter- 
mediate energy) results in a merged system. The global 
observables are plotted in Fig. [2j It may not be obvious 
at first why the fragments should split for lower ener- 
gies and merge for higher ones. But the estimate for the 
occupied phase-space volume Vk(t) presented in subplot 
(c) indicates that V increases with energy. Therefore the 
fragments' average distance in phase space is larger, while 
in compensation they can come closer to each other in co- 
ordinate space. However, this behavior is also dependent 
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Video 1| (color online) Two-dimensional z-fc z -slice from 
the full six-dimensional Wigner distribution /^(r, k) for a 
central 40 Ca+ 40 Ca collision with a center-of-mass energy of 
£c.m. = 160 MeV. 
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|Video 2| (color online) Same as Video [5] with a center-of- 
mass energy of E c . m . = 240 MeV. 



on the particular Skyrme force used and the presence of 
time-odd terms, which is discussed in the next subsec- 
tion. 



C. Influence of time-odd terms 

Skyrme energy-density functionals are calibrated to 
ground state properties of even-even nuclei [PEU2T)] . This 
leaves the choice of the time-odd terms in the functional 
(current j 2 , spin-density s 2 , spin kinetic energy density 

-H- 

T, and the spin-current pseudotensor J) largely unspec- 
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FIG. 2. Global observables /3*(t) (a), j k (t) (b), V h (t) (c), 
and Qf{t) (d) are shown for a central 40 Ca+ 40 Ca collision 
with fixed Skyrme interaction SkI4. Each curve corresponds 
to a different center-of-mass energy. 



ified [23]. Galilean invariance requires at least some of 
these terms to be present depending on the presence of 
the associated time-even term, e.g. j 2 for the (pr — j 2 ) 
combination. In our calculations we always include the 
time-odd part of the spin-orbit interaction. In order to in- 
vestigate the effects of the remaining time-odd terms, we 
have compared different choices by using a single Skyrme 
parametrization and the same test case. We choose the 
force SLy4 and start with the minimum number of time- 
odd terms which is needed to ensure Galilean invari- 
ance [S3]. In the next stage, we include also the spin- 
density terms proportional to s 2 . Finally, we also add the 
combination which includes the tensor spin-current term 
(s • T — J 2 ) . A comprehensive notation of the full Skyrme 
functional can be found, e.g. in Ref. [26]. As shown in 
Fig. [3| at least for the quantities j3 k and Q 20 , varying 
these time-odd terms has a very little effect in the ini- 
tial contact phase and the dynamical behavior becomes 
somewhat different only in later stages of the collision. 
On the other hand, small differences near the threshold 
energy (the highest collision energy for a head-on colli- 
sion that results in a composite system. At higher en- 
ergies the nuclei go through each other) can have large 
long-term effects on the outcome of the collision. For 
example a small difference in dissipation may be enough 
to influence the decision between re-seperation or form- 
ing a composite system. We have also checked a broader 
range of collision energies from the fusion regime up to 
deep inelastic collisions. The interesting quantity is the 
loss of fragment kinetic energy between the entrance and 
exit channels. It was found that the spin terms contribute 
small changes to this loss which can go in both directions, 
less dissipation near fusion threshold and more dissipa- 
tion above. Subplot (b) of Fig. [3] shows the effects near 
the Coulomb barrier where spin terms reduce dissipation. 



IV. SUMMARY 

Wc have presented a geometrically unrestricted frame- 
work to study nuclear dynamics within TDHF in the full 
six-dimensional phase space. The impact of different re- 
action parameters on the outcome of a heavy-ion collision 
was studied in detail for the 40 Ca+ 40 Ca system. We find 
that the occurrence of transparency is clearly reflected in 
the global asymmetry of the Wigner momentum distri- 
bution. The surprising result that in some cases the sys- 
tem merges at higher energies and shows transparency 
at lower ones can be related to the interplay between 
momentum- and configuration-space volumes which is a 
reflection of the Pauli principle. It is also interesting that 
the two distributions in phase-space never truly com- 
bine to form a single distribution. This clearly indicates 
that two-body collisions will be necessary to achieve true 
equilibrium as the reaction proceeds to longer contact 
times. The detailed degree of relaxation found depends 
on energy and also the properties of the Skyrme force, 
where especially the effective mass seems to be impor- 
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tant. The presence of additional time-odd terms in the 
Skyrme functional appears to have a complex impact on 
the outcome of a collision as well. In this paper only one 
non-central collision was studied. A systematic investiga- 
tion of impact parameter and energy dependence (fusion, 
deep-inelastic reactions) as well as even heavier systems 
would be highly interesting but is beyond computational 
feasibility at the moment. 
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FIG. 3. Global observables /3 k (i) (a), and Q 2 r °(t) (d) are 
shown for a central 40 Ca+ 40 Ca collision with fixed Skyrme 
interaction SLy4. 
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